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Abstract 

The creation of electron-positron pair by a photon and the bremsstrahlung 
of an electron in a medium are considered in high-energy region, where in- 
fluence of the multiple scattering on the processes (the Landau-Pomeranchuk- 
Migdal (LPM) effect) becomes essential. The pair photoproduction proba- 
bility is calculated with an accuracy up to the " next to leading logarithm" . 
The integral characteristics: the radiation length and the total probability 
of pair photoproduction are analyzed under influence of the LPM effect, 
and the asymptotic expansions of these characteristics are derived. 
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1 Introduction 



When a charged particle is moving in a medium it scatters on atoms. With 
probability ~ a this scattering is accompanied by a radiation. At high energy 
the radiation process occurs over a rather long distance, known as the formation 
length l c : 

lc = 1 ° 2 Q2' ^° = — 2~ ' (■'■•-'■) 
1 + m z oo 

where uo is the energy of emitted photon, e(m) is the energy (the mass) of a 
particle, e' = e — oo, d c is the characteristic angle of photon emission, the system 
h — c — 1 is used. The spectral distribution of the radiation probability per 
unit time inside the thick target (the boundary effects are neglected) can be 
obtained from the general formula for the spectral probability derived in the 
framework of the operator quasiclassical method (see Eqs.(4.2)-(4.8) in 0). It 
can be estimated as 

a A 2 (/ c ) duo 
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where a = e 2 = 1/137, A 2 (Z C ) is the mean square of momentum transfer to a 
projectile from a medium (or an external field) on the formation length l c . 

If the angle of multiple scattering on the formation length $ s = \f^t c is small 
comparing with the angle I/7 (7 = e/m is the Lorentz factor), then one can 
consider scattering as a perturbation and perform the decomposition over "the 
potential" of a medium. The radiation probability in this case is the incoherent 
sum of the radiation probabilities on isolated atoms of a medium defined by 
the Bethe-Heitler formula. One get from (|1.2|) for the spectral probability of 



radiation per unit time in the case $ s I/7 ($ c = 1/7, A 2 
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At an ultrahigh energy it is possible that ?9 S ^> 1 /j. In this case the character- 
istic radiation angle (giving the main contribution into the spectral probability) 
is defined by the angle of multiple scattering i} s . The self-consistency condition 
is 



' c ^ . 2 • 



From the condition (|1.4j) we find 
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In this case one get from (|1.2j ) for the estimate of the spectral radiation proba 
bility per unit time 
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where Z is the charge of the nucleus, n a is the number density of atoms in the 
medium, A c = 1/m = (h/mc) is the electron Compton wavelength, a s is the 
screening radius of the atom. 

So, the formula (|1.2|) gives the general description of the radiation process 
in terms of the mean momentum transfer valid both in a medium and in an 
external field, while formulas ( |1.3|) and ( |1.6|) describe the process probability in 
the particular regimes in a medium. 

Landau and Pomeranchuk were the first who showed that if the formation 
length of bremsstrahlung becomes comparable to the distance over which the 
multiple scattering becomes important, the bremsstrahlung will be suppressed 
Migdal || developed the quantitative theory of this phenomenon. 

New activity with the theory of the LPM effect (see [ffl, |J, ||) is connected 
with a very successful series of experiments performed at SLAC recently (see @ , 
M). In these experiments the cross section of the bremsstrahlung of soft photons 
with energy from 200 keV to 500 MeV from electrons with energy 8 GeV and 
25 GeV is measured with an accuracy of the order of a few percent. Both LPM 
and dielectric suppression are observed and investigated. These experiments 
were the challenge for the theory since in all the mentioned papers calculations 
are performed to logarithmic accuracy which is not enough for description of 
the new experiment. The contribution of the Coulomb corrections (at least for 
heavy elements) is larger then experimental errors and these corrections should 
be taken into account. 

Authors developed the new approach to the theory of the Landau- Pomeranchuk- 
Migdal (LPM) effect @ in which the cross section of the bremsstrahlung process 
in the photon energies region where the influence of the LPM is very strong was 
calculated with a term oc 1/L , where L is characteristic logarithm of the prob- 
lem, and with the Coulomb corrections taken into account. In the photon energy 
region, where the LPM effect is "turned off", the obtained cross section gives the 
exact Bethe-Heitler cross section (within power accuracy) with the Coulomb cor- 
rections. This important feature was absent in the previous calculations. Some 
important features of the LPM effect were considered also in [|K| , |lTJ , Jl2| . 

The crossing process for the bremsstrahlung is the pair creation by a photon. 
The created particles undergo here the multiple scattering. It should be empha- 
size that for the bremsstrahlung the formation length ( |1.4| ) increases strongly 
if uj <C e. Just because of this the LPM effect was investigated at SLAC at a 

relatively low energy. For the pair creation the formation length L = — — - 

m 2 uj 

attains maximum at e = uj/2 and this maximum is l Pi max — (uj/2m)X c . Be- 
cause of this even for heavy elements the effect of multiple scattering becomes 
noticeable starting from u ~ 10 TeV. Nevertheless it is evident that one have 
to take into account the influence of a medium on the pair creation and on the 
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bremsstrahlung hard part of the spectrum in electromagnetic showers being cre- 
ated by the cosmic ray particles of the ultrahigh energies. These effects can be 
quite significant in the electromagnetic calorimeters operating in the detectors 
on the colliders in TeV range. 

In the present paper both the spectral probability and the integral proba- 
bility of the pair creation are calculated within an accuracy up to "the next 
to logarithm" and with the Coulomb correction taken into account (Sec. 2). In 
Sec. 3 the radiation length is calculated under influence of the LPM effect. The 
total probability of photon radiation is considered also. In the Appendixes the 
technical details of calculation are given. 

2 Influence of multiple scattering on pair cre- 
ation process 

The probability of the pair creation by a photon can be obtained from the prob- 
ability of the bremsstrahlung with help of the substitution law: 

u 2 du> — > e 2 de, lo — > — 10, e — > — e, (2.1) 

where u is the photon energy, e is the energy of the particle. Making this 
substitution in Eq.(2.12) of || we obtain the spectral distribution of the pair 
creation probability (over the energy of the created electron) 

dW p _ 2auf lm ^ ^ _ + ^ ^ _ ^ ^ (2 2) 
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G = H+l, H = p* -iV{g), 

2nZ 2 a 2 ee'n a r , a 2 2 
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where C = 0.577216 ... is Euler's constant, n a is the number density of atoms 
in the medium, q is the coordinate in the two-dimensional space measured in 
the Compton wavelength A c , which is conjugate to the space of the transverse 
momentum transfers measured in the electron mass m. The mean value in (|2.2|) 
is taken over the states with the definite value of the operator q (see [[J], Sec. 2). 
The contribution of scattering of the created electron and positron on the atomic 
electrons can be incorporated into the effective potential V^(^) by substitution 

1.88 + f(Za) 



Q -> Qef = <2(1 + -^r), a S 2 — > a ef = a s2 exp 



(2.4) 
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The potential V(g) in Eg. ( |2.3|) we write in the form 

V(g) = V e (e)+v( e ), V c (g) = gg 2 , q = QL C , 

Wft)=ln j|, v{e)= .^^ +2C y ( 2 .5, 

where the parameter g c is defined by the set of equations: 

g c = l for 4QL 1 <1; 4Q^ln^ 2 - = l for 4QL X > 1, (2.6) 

where L\ is defined in Eq.( |2.3| ). The parameter q c ~ l/p c is determined by 
the characteristic angles of created particles with respect to the initial photon 
momentum (or the corresponding momentum transfers). In accordance with 
such division of the potential we present the propagators in the expression (|2.2|) 
as 

G — Gq = G — G c + G c — Gq (2-7) 

where 

G c = H c + 1, G = H c + 1 - iv, H c = p 2 - iqg 2 

This representation of the propagator G^ 1 permits one to expand it over the 

"perturbation" v. Indeed, with an increase of q the relative value of the per- 

v 1 

turbation diminishes (— ~ — ) since the effective impact parameter diminishes 

► c L c 

and, correspondingly, the value of logarithm L c in (|2 . 5|) increases. The maximal 
value of L c is determined by a size of a nucleus R n 

a 2 

L mm =\*-±~2L X , (2.8) 

where a s = a S 2exp(/ — 1/2) = lll.Z , ~ 1 / 3 A c . When g c R n one cannot consider 
the potential of a nucleus as the potential of a point charge. In this case the 
expression for the potential V(g) has been obtained in [|], Appendix B 

V{q) =qQ 2 {L max - 0.0407). 

The matrix elements of the operator G^ 1 was calculated explicitly in 0: 

oo 

< ^i I CfT 1 1 ^2 >= i J dte~ lt < Q 1 \exp(—iH c t)\Q2 >, 

o 

< 0i | exp(-zH c t)|£ 2 >= K c (g 1 ,Q 2 ,t) 



v r w 



exp < — 

Ani sinh ut 14 



2 

\q\ + q\) coth vt - . Q 1 Q 2 
smh ut 



(2.9) 
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where v = 2yiq. Substituting this expression in the formula for the spectral 
distribution of the pair creation probability (|2.2|) we have 



dW c v 
de 



Im ®Jv) 



dte 



si 
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V sinh z 



sinh z 
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where z = vt, p 



i/(2v), ip(x) is the logarithmic derivative of the gamma 
function. Some details of the derivation of the last line can be found in Appendix 
A (see (|A.1|) - ([A.8|) ). This formula gives the spectral distribution of the pair 
creation probability in the logarithmic approximation which was used also by 
Migdal || . It should be noted that the parameter g c entering into the parameter 
v (see Eqs.(|2.3D and (|2.5|)) is defined up to the factor ~ 1, what is inherent 
in the logarithmic approximation. However, below we will calculate the next 
term of the decomposition over v (q) (an accuracy up to the "next to leading 
logarithm") and this permits to obtain the result which is independent of the 
parameter g c . It will be shown that the definition of the parameter g c in Eq.(2.6) 
minimizes corrections to ( |2.10| ) practically for all values of the parameter g c . It 
should be emphasized also that here the Coulomb corrections are included into 
the parameter v in contrast to ||. 

Let us expand the expression G^ 1 — G^ 1 over powers of v 



G' 



G r 



G~ {iv)G~ + G~ {iv)G~ {%v)G~ + 
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Substituting this expansion in ( [2.6|) and then in (|2.2| ) we obtain the decompo- 
sition of the probability of the pair creation. Let us note that for q 1 the 

dw; /R . 

sum of the probability of the pair creation J _ ( p.l0|) and the first term of the 



expansion ( p. 11 ) gives the Bethe- Heftier spectrum of electron of created pair, see 
below (|2.22j ). At q > 1 the expansion (|2.11| ) is the series over powers of — . It 



is important that the variation of the parameter g c by a factor order of 1 has an 
influence on the dropped terms in (|2.11|) only. 

In accordance with ( |2.7p and ( |2.11[ ) we present the probability of radiation 
in the form 

dW n dWZ dWl dWl 



de 

The probability of pair creation 
with allowance for 



de 
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de 



de de 
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+i J dh J dt 2 J dt 3 e- l{h+t2+h) J d 2 Ql J d 2 g 2 K c (0, q x M)v{q 1 )K c {q 1 , g 2 ,t 2 



xv(g 2 )K c {g 2 ,0,t 3 ) + 



where the matrix element K c is defined by 



dWl 

The term — — — 

de 



(2.13) 
in (£3) 



corresponds to the first term (linear in v) in fl2.13| ). Substituting (|2.9| ) we have 



dWl 
de 
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Substituting in ( 2.14 ) the explicit expression for v(g) and integrating over d 2 g 
and <i(ii — 1 2 ) we obtain the following formula for the first correction to the pair 
creation probability 



de 
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g(z) = z cosh z — sinh z, t = t\ + t 2 , z = vt 
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(2.15) 



here Li 2 (x) is the Euler dilogarithm. Use of the last representation of function 
G(z) simplifies the numerical calculation. 

As it was said above (see (|2.6p, (|2.8|)), g c = 1 at 



v\ = AQL X < 1 (g = QLi 



(2.16) 



If the parameter \v\ > 1, the value of g c is defined from the equation (|2.6|) . Then 
one has 



In g\ + In - = ~ ln(^4QL c ) - = -z*, ^4QL C = 1. 



z 2 _vcc ""' _w "4 4' cc ---" - ^ 2 - 17 ^ 
So, we have that the factor at g(z) in the expression for fi(z) in ( |2.15 ) can be 
written in the form 



v 



(In gi + In In sinh z - C) -> (In - + lnfiWW^ - 1) 
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where 
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L(Qc), 



(2.19) 



is the Heaviside step function. So, we have two representation of \u\ de- 



Vq. 



pending on g c : at g c = 1 it is |z/| = z/i and at C < 1 it is 

When the scattering of created particles is weak (v\ <C 1), the main contri- 
bution in Q2.15D gives the region where z <C 1. Then 



F P {v) 



3 v 3 



C-ln(it)), 



si 



L -> Li. 
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Substituting the expansion (|C.1| ) into Eq. (|2.10| ) we find the corresponding asymp- 
totic decomposition of the function Q p (v) 



$ p (z/) ~ Si 



if + yr 



31z/ 6 
126 



+ s 2 



7/ 

y 



2z/ 4 16z/ 6> 
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Combining the results obtained ( 2.20| ) and ( 2.21 ) we obtain the spectral distri- 
bution of the pair creation probability in the case when the scattering is weak 

(M«D 
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where L\ is defined in 
28Z 2 a 3 n 



where 



9m 2 



. Integrating Q2.22 ) over e we obtain 

3312 uj 2 \ 
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(2.24) 



Note that in gold u Q = 10.5 TeV. This is just the value of photon energy starting 
with the LPM effect becomes essential for the pair creation process in heavy 
elements. If one omits here the terms oc v\ and oc {uj/ujq) 2 these expressions 
coincide with the known Bethe-Heitler formula for the probability of pair creation 
by a high-energy photon in the case of complete screening (if one neglects the 
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contribution of atomic electrons) written down within power accuracy (omitted 



m 



terms are of the order of powers of — ) with the Coulomb corrections, see e.g. 



Eqs.(19.4) and (19.17) in P 



UJ 



The pair creation spectral probability dW/dx vs x = e/uj is shown in Fig.l 
for different energies. It is seen that for uj = 2.5 TeV which below ujq the 
difference with the Bethe-Heitler probability is rather small. When uj > ujq 
there is significant difference with the Bethe-Heitler spectrum increasing with uj 
growth. In Fig.l are shown the curves (thin lines 2,3,4) obtained in logarithmic 
approximation dW^/de ( p.lOj ), the first correction to the spectral probability 
dWp/de ( |2.15|) , curves c2,c3,c4 and the sum of these two contributions: curves 
T1,T2,T3,T4. It should be noted that for our definition of the parameter q c 
( |2.6| ) the corrections are not exceed 6% of the main term. The corrections are 
maximal for ~ 3. 

The total probability of pair creation in the logarithmic approximation can 
be presented as (see fl2.10|) ) 



wi_ 



_9_^o Im f dy 
14 u m J 2/(1-2/) 



In p — ip [ P + - 



1 -2y + 2y 2 ) lip(p) -lnp + 



2 P/ 
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where 



bs 



y/y^-y)' 




is the Bethe-Heitler probability of pair photoproduction in the logarithmic 
approximation. The total probability of pair creation W£ in gold is given in Fig. 2 
(curve 2), it reduced by 10% at uj ~ 9 TeV and it cuts in half at uj ~ 130 TeV. 

At vq 1 the main term of the function F p (u) (see (|2.15|) and ( 2.19|) ) can be 
written in the form 



dz 



sinh z 



[s 1 f 1 (z)-2is 2 f 2 {z)}. 



(2.26) 



Integrating over z we obtain 



7T 



-Im FJv) = ~(si 



S2) + ^( ln2 



7T 



(2.27) 



where we take into account the next terms of the decomposition in the term 
cx s 2 - Under the same conditions (u ^> 1) the function Im 3> p (f) ( p.10 ) is 



Im %(u) = -(si - s 2 ) + -J=s 2 . 



[2.28) 
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dW} 

Thus, at Pq ^> 1 the relative contribution of the first correction — is defined 

de 

by 



dW[ 
~dWj 



2L r 



In 2 



c + i 



0.451 



(2.29) 



In this expression the value r with the accuracy up to terms ~ 1/1% doesn't 
depend on the energy:L c ~ L\ + \n{u /ojq) /2. Hence we can find the correction to 
the total probability at uo uq. The maximal value of the correction is attained 
at uj ~ lOc^o, it is ~ 6% for heavy elements. 

When the parameter v\\ is not very large (u < 10 3 , g c > R n , see Q2.8|) ) one can 
solve the equation v\\q\ = \ (|2.17|) using the method of successive approximations. 
In the first approximation we have 



1 

cosh£ 



L ~ 



U 1 1 + 



1 



Li 
ln^ 



2 In cosh £ 



(2.30) 



It should be noted that the relative error in the expression for L c at q c > R n is 
less than In 2/(41^) < 2.5%. Here we introduce variable £ 



£ 1 /, i . ££' 

- = -(1 + tanhf), 



1 



UJ 



4cosh^ 



-oo < £ < oo. 



(2.31) 



Substituting the terms oc v§ in the asymptotic formulas ( 2.28 ) and ( |2.27| ) into 
the first line of Eq. ([2.22|) we obtain expression which contain the integral of the 
type 



cosh£ 



2coshX 



37T 

[A + B In cosh 0] = — 



A + B (ln2 + - 



. (2.32) 



Using this result we obtain the total probability of pair creation under strong 
influence of multiple scattering (vq j$> 1, but not very large) 
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4v / 2 \/uJUJo 
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2C 



1 

3 



4Li 
0.08303 - 



In — + D 

UJ 
1 



12.04 



(2.33) 



It should be noted that only the main term of the decomposition (oc u ) can 
be used in Eqs. fl2.27l) and (|2.28|) for the calculation of the total probability of 
pair creation. In the interval uj ujq the contribution into the correction terms 
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gives also the region where cosh 2 £ ~ lo/loq, where the parameter is ~ 1 and the 
expansion used in Eq. ( |2.26| ) is ineligible. The next terms (without corrections 
oc l/-^i) are found in Appendix A (Eq.( |A.12] )), so we have 



3a m 2 
4 a/2 s/loloq 



1- Y(41n2-1) 



,co vr^ + 1 ,„ +D 
to 18 to ALi V oo 



[2.34) 



In terms of the Bethe-Heitler total probability of pair creation this result is 



0.836, 



0.548^+ ' 



to 



„ . In — + 0.274 

to ALi \ loq 



(2.35) 



3 Influence of the multiple scattering on the 
bremsstrahlung 



The spectral radiation intensity obtained in U (see Eq.(2.39)) has the form 



where 



dl = codW = ^^-Im 
27r(l — x) 



$(u) - 



2L, 



-F{v) 



LO 
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(3.1) 
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( 1 



V sinh z z 
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IVT 2 



1 



sinh z z 2 
1 \ 



ri ( Inp - i) ( p + J + r 2 hp(p) - hip + . , 



F{v) 



-Arr- [rih(z) - 2ir 2 f 2 (z)} 
smh z 



z 

j 

V 



T\ — X 



r 2 — 1 + (1 — x) 



(3.2) 



where z = vt, p = %j (2v\ ip(x) is the logarithmic derivative of the gamma func- 
tion. Some details of the derivation of the second line can be found in Appendix 
A (see (pOP- ([OP ). The functions f x (z) and f 2 (z) are defined by Eq.(frjg), 



v =w/ , V Q 



e 1 — x 



e x 



e = m [&TxZ 2 ot 2 n a \ i c Li 



L c ~ L x ( 1 + " 1) 



(3.3) 



Note, that the parameter e is four times smaller than the parameter loq defined 
in Eq. (|2.24D. The LPM effect manifests itself when 



v>i(x c ) 



E + E 



(3.4) 
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The formulas derived in |J and written down above are valid for any energy. 
In Fig. 3 the spectral radiation intensity in gold (sq = 2.5 TeV) is shown for 
different energies of the initial electron. In the case when e <C 6q (e = 25 GeV 
and e = 250 GeV) the LPM suppression is seen in the soft part of the spectrum 
only for x < x c ~ e/eq <ti 1 while in the region e > Eq (e = 2.5 TeV and 
e = 25 TeV) where x c ~ 1 the LPM effect is significant for any x. For relatively 
low energies e = 25 GeV and e = 8 GeV used in famous SLAC experiment 0, 
[§] we have analyzed the soft part of spectrum, including all the accompanying 
effects: the boundary photon emission, the multiphoton radiation and influence 
of the polarization of the medium. The perfect agreement of the theory and 
data was achieved in the whole interval of measured photon energies (200 keV< 
uj <500 MeV), see the corresponding figures in P|,[p!0|,[p|. It should be pointed 
out that both the correction term with F(v) and the Coulomb corrections have 
to be taken into account for this agreement. 

In the case e <C Eq in the hard part of spectrum (1 > x ^> x c ) the parameter 
v\ ~ x c /x <C 1 and the contribution into the integral (|3.2|) give the region z < 1. 
Using the decomposition flC.lD we find (compare with (|2.21| ), (|2.22j) ) 



Im $(z/) ^ 
-Im F(u) 



v. 



rv 



31 
21 



r 2 - 



1 - 



16 
21 



^i 4 I 



9 



(^-r^^a + OK 4 )). 



(3.5) 



In the last formula, which presents corrections ~ 1/Li we restricted ourselves to 
the main terms of expansion. Substituting into (|3.1|) we have 



dl _ 2Z 2 a 3 n a e 
dx 




+2r 2 



3m 2 



+ 



6 



x 



(3.6) 



Note that if neglect here the terms cx x\jx 2 we obtain the Bethe-Heitler intensity 
spectrum with the Coulomb corrections. 

In the case e > 6q the intensity spectrum differs from the Bethe-Heitler one 
at x ~ 1 also. When e ^> e one can use the asymptotic expansions ( |2.27| ) and 
([2.28 ) in the interval not very close to the end of the spectrum (x — 1): 




(3.7) 
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Now we turn to the integral characteristics of radiation. The total intensity 
of radiation in the logarithmic approximation can be presented as (see Q3.1|)) 



i 



+ 



dx 



x 



g V 1 - x 



(2(1 - x) + x 2 ) 



Ul>{p + 1) - 4> (p + ^)) + 2 J xdx(ip(p + l) 



— hip) 



(3.8) 



where 



9V 

P= Y> 



V 



x 



x 



9 



exp \i- 




L° rad is the radiation length in the logarithmic approximation. The relative energy 
losses of electron per unit time in terms of the Bethe-Heitler radiation length L° rad : 

—L° rad in gold is given in Fig.2 (curve 1), it reduces by 10% (15% and 25%) at 

e ~ 700 GeV (e ~ 1.4 TeV and e ~ 3.8 TeV) respectively, and it cuts in half 
at uj ~ 26 TeV. This increase of effective radiation length can be important in 
electromagnetic calorimeters operating in detectors on colliders in TeV range. 
The contribution of the correction terms was discussed after ( |2.29| ). It is valid 
for the radiation process also. 

The spectral distribution of bremsstrahlung intensity and the spectral dis- 
tribution over energy of created electron (positron) as well as the reduction of 
energy loss and the photon conversion cross section was calculated by Klein |T^| . 
[§] using the Migdal formulas. As was explained above (after Eq. ( |2.10|) ) we 
use more accurate procedure of fine tuning and because of this our calculation 
in logarithmic approximation differs from Migdal one. We calculated also the 
correction term and include the Coulomb corrections. For this reason the results 
shown here in Figs. 1-3 are more precise than given in fI5| , ||. 

In Eqs. (|3.6j ) and (|3.7p we can use the main terms of decomposition only. The 
main term in ( |3.6| ) gives after the integration over x the standard expression for 
the radiation length L ra d without influence of multiple scattering. The correction 
term is calculated in Appendix C (see ( |C.9| )) where we need to put |/3| 2 = eq/e ^> 
1 



am 



1 



1 



Atxeq V ' 9Li 
2Z 2 a 3 n a L 1 



J rad 



4ir £ 
15 e 



r-i 

^rad 



1 - 



47T e 

15 Eq 



9Li 



(3.9) 



The integration over x of the main term in (|3.7| ) gives (terms oc y£o/e in the 
square brackets are neglected) 



9irZ 2 a 3 n av fEE^ T 
to , 7= n Li 



Ay/2- 



1 + 



4Li 



. e 46\ 
ln ^"27j +r ° 
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T 



2Li 



7T 



In2-C + -]. 



(3.10) 



The corrections (without terms oc 1/Li) to ( |3.10| ) are calculated in Appendix B 
(see Eq. QB.lip ). The complete result is 



9am [~e 



32y/2V e 



9 v J \ e 54 e 



1 



4Li 
I 



e 7r 46 

In— + 21n2-2C + -- — 
e 2 27 



eL r , 



rad 



5 /£q 

2V e 



1 - 2.37 



4.57^+ ' 



4Li 



In 0.3455 



(3.11) 



Although the coefficients in the last expression are rather large at two first terms 



of the decomposition over yEo/s this formula has the accuracy of the order of 
10% at e ~ lOeo- The integral probability of radiation was calculated in [Til 



W 



3L r 



ad 



ln^ + C 2 

e 



Co = 2C 



+ \nz 



1 cosh z ' 
z 3 sinh 3 z 



1.96 (3.12) 



In the case e ^> eo we can calculate the integral probability of radiation start- 
ing with Eq. ( |3.7D . Conserving the main term, dividing it by xe and integrating 
over x we find 



_ llnZ 2 a 3 n a [ej 
w o — — — - — \l — L i 



2V2: 



1 /, £ 8 

1 H In 1 

4Li V e 11 



(3.13) 



The correction terms to Eq. (|3.12|) are calculated in Appendix B (see Eq.( B.13| )). 
Substituting we have 



W 



Ham 2 



16v^i^ 



1-^(21*2 + 1) 

11 V e 6 e 



— (in— + 21n2-2C + - 
4L X V £n 2 11 



llnZ 2 a 3 n a [eq 
2^2m 2 



I - 1.23 J— + 1.645— + ' 



4Lj 



In — + 2.53 



(3.14) 



Ratio of the main terms of Eqs. (|3.11|) and (|3.14[) gives the mean energy of 
radiated photon 

u = —e ~ 0.409^. (3.15) 
22 
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A Appendix 

We consider the integral which represent the integral probability of pair photo- 
production (see Eqs. fl2.10D and ( |2.31| )) 



oo oo 



n(a) = J cf£ J dzexp(— a z cosh 







1 



1 1 



sinh z z 



a cosh £ y 2 cosh £ / V sinh z z 



(A.l) 



Integrating by parts (over z) the second term of the integrand in ( A.l|) we have 



Ma) 



)oo oo 
+ J d£ J dzexp(—az cosh£) 
n n 



1 



+1 — coth ,2 



a J cosh £ \ 2 cosh 
o s v 

1 







1 



1 — coth z H — 

z 



2 coslT £ 

The functions entering in (|A.2p we present as 
1 



sinh z 
(A.2) 



2j]exp(-(2A;- l)z), coth z- 1 = 2^exp(-2fcz) 
fc=i fc=i 



sinh z 

Let us consider the integral entering ( A.2| ) 

d£ 



(A.3) 



7riiaj 



-2-F!(acosh£), 

cosh £ 



Fi(acosh£) = ^ dzexp(— az cosh£) ^1 — coth^ H — ^ (A. 4) 

o 

To avoid a divergence of the individual terms in the integral over z we put the 
lower limit of the integration 5 — > 0. Using (|A.3|) we obtain 



Fi (x) = lim 



Using the expansions 



~ exp(-2^) 
-Ei(-Jx) - ^ - + ^ 



x 



k=l 



\ k(2k + x) 



(A.5) 



-Ei(-fe) = -ln(«5z)-C, - f] 6Xp( 2M) = ln(l - exp(-2£)) = In 2(5, 



fc=i 



=^ - + 1 -In 



.r 



(A.6) 
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where ip(x) is the logarithmic derivative of the gamma function, and taking 
integrals over £ we have 



4 a 00 1 n 00 
7ri(a) = In — C-H — V TT7-a V — . 

The formula ( |A.2| ) contains also the integral 



00 

^2(0) = y d^F 2 (a cosh£) 


00 1 

-FU^) = / dzexp(-zx) ( 1 — cothz H 

J V smh 2; 



1 , yj2k + a + V%k 



In 



2a 



(A.7) 



^ - + 1 -V 



x + 1 



(A.8) 



Transposing the integration order and using (|A.3|) we find 

00 y 

n 2 (a) =2J2 dzK (az) (exp(-(2A; - l)z) - exp(-2A;z)) , (A.9) 
k=i{ 

where K (x) is the modified Bessel function. Taking here integrals we have 



?r 2 (a) = 2£ 



k=l 



^j(2k - l) 2 - a : 



dn 



2Jfe - 1 + V( 2A; - !) 2 - fl2 



In 



2fc + - a 2 



V4A; 2 - a 2 a 
Substituting flOD and flA.lOQ into Eq. flA^) we obtain 



(A.10) 



„, > 37r 1 /. a 

n ( a ) = -r+ol ln - + 1 + c 

8a 2 



7r 3 a 
l8~ 



fe=i 



2fcV4fc 2 - a 2 



In 



V2A; + a + 



2a 



(2Jfe - l) 2 - a 2 



In 



2 A; - 1 + J (2k - l) 2 - a 2 



2 , 2k + V4A; 2 - a 2 
in 



2 - a 2 



a 



(A.ll) 



This expression is particularly convenient at \a\ < 1. In the case \a\ <C 1 the 
first three terms of the decomposition are (a = |a| exp(?7r/4)) 



i7T 7T 

Im 11(a) ~ — ;=— + - (1 -41n2) 



8^^ 8 



7T 3 |a| 



48^ 



(A.12) 



16 



B Appendix 

Here we consider the asymptotic behavior of the radiation integral characteristics. 
The integral intensity (the radiation length) can be presented as (see ( |3.1|) and 

my. 



1 oo 

./. (./ </. 



1 — X 



+ —(l + (l-x) 

Pv 



dzexp(—(3rjz) 

(-2 - 

Vsinh 2 



x 



( 1 



1 



Vsinh 



z z 



z z z 



X 



1 — X 



(B.l) 



Integrating by parts (over z) the second term of the integrand (oc l/((3r])) in 
(|B.1| ) we have 



P((3) = ~J ^^(l + il-x) 2 ) dx + JJ^LJdzeKpi-foz) 



X 



~) + (l + (1 - X) 2 ) (l - COtllZ + -) 

sinnz z) v J V z) 



9tt 
16/3 



where 



+ P 1 (P)+P 2 (P), 



(B.2) 



r x dx r f I 

= / / dzexp(-(3r]z) 1 - coth z H 

J 1 — x J \ sinh z 

o o 

1 oo i 

^2(/3) = 2 J %dx J dzexp(—f3riz) ^1 — cothz H — 



o o 

We split the function P{(3) into two functions: 

P 1 (P)=P n {/3) + P 12 (l3), 



(B.3) 



P n (/5) = / — — [ dzexp(-/3riz) (l - coth z H ^ 

J 1 — x J V smh z 

o o 

1 oo ^ 

P 12 (/3) = - y "(1 + x + x 2 )dx J dzexp(-f3r]z) M -cothz+ — 



sinh z 
o o 

Transposing the integration order in Pll(f3) we get the integral over x 



(B.4) 



o o r a 

oo 

= -2In(/3z) + /" ln(/?V + y 2 ) exp(-y)dy 



(B.5) 
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In the limit <C 1 we have discarding terms ~ (5 2 

oo 

P xx {0) = -2 J dz(\n((3z) +C)(l- cothz + — — 
o 

We use Eq.( |A.3| ) in the calculation of the integral over z in Pi2(/3) 

7 / i \ 

/ dzexp(-priz) 1 - coth^ + — — = 2 V 
J \ sinhz/ ZT 1 , 



(B.6) 



o 



k= \\2k - 1 + p V 2k + (3r]J 

(B.7) 



21n2-^C(2). 

Taking into account that j3 = \/3\ exp(i7r/4) we have for ImP(/3) at \/3\ <C 1 

1 



ImPn(/3) = -| J dz (l - coth^ + - 



sinh z 



ImP 12 (/?) 



v/2 



1 — X 



(1 + x + x 2 )dx 



-vrln2 
19tt 



|/?|C(2). (B.8) 



We will use Eqs.( A.4 )-( A.6 ) in the calculation of the integral over z in P2{/3) 
Eq-O 



7 / 1\ 2 00 
/ ofo exp(— /3r]z) 1 — coth — I = In — C + 



£77 



ln|--C+W(2) 



(B.9) 



As a result we get 



ImP 



1 

m = 2 / 



-H^/^H-i + ^ c(2) (Bio) 



Substituting (p.8p and ( |B.10| ) into ( |B.2| ) we obtain 
ImP(/3) = 



^ - 7 rln2-^ + ^C(2)+0(/? 2 ) 



16V2|/3| 



4 l6y/2 



(B.ll) 



In the calculation of the total probability of radiation one have to make the 
substitution in ( p.l| ) 



xdx 
1 — x 



dx 
1 — x 
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Then 



(B.12) 



where 



Im t x {(5) = -7rln2 + ^=|j9|C(2) / ^7^(1 + = - nln2 + £^ I/3|C(2) 



8^ 



Im t 2 (/3) 
Im 



7T 1 



2 2v/2 

117T 7T 



X 



1 — X 

lbr 



-da; 



7T 7T 



2 2x/2 



l^iC(2), 



WW 2 -l^ + ^IC(2) + 0(^) 



(B.13) 



C Appendix 



Here we consider the asymptotic behavior in the region where the LPM effect is 
weak. In this region in Eq.(A.l) the parameter \a\ ^> 1 and in the integral over 
z the interval z -C 1 contributes, and we can use expansions 

1 1 „ 7z 3 31z 5 



sinh z 
1 



6 + 360 



1 

3 



z 
15 



15120 

2z 4 



+ . . . 



sinh'z z 2 3 15 189 

Substituting these expansions into (|A.1|) and taking into account that a 
\a\ exp(i7r/4) we get 



(C.l) 



Im 11(a) 



d£ 



1 



31 



6|a| 2 cosh £ 126|a| 6 cosh £ 



+ 1- 



1 



2coshX 



16 



7 



184 



X ^3|a| 2 cosh 2 £ m\a\^co^i} 1 °'~ 12 ' 1 * AO '~ 14 1 " ( °' 2) 

We turn to Im T(J3) Eq.( gg) at > = \/3\ exp(i7r/4)). The integral 
over z coincides with this integral in ( |C.1| ). The integral over x gives for x ~ 1 the 
same structure as in ( p.2| ): the main term ~ 1/|/?| 2 and the correction ~ 1/|/3| 6 . 
In the region x ~ 1/|/5| 2 where the influence of the LPM effect is significant, 
the correction is proportional to the phase space: J xdx ~ 1/|/3| 2 . This is the 

main correction. Because of this we split the integration interval over x into two 
intervals: 1) < x < Xq and 2) Xq < x < 1, where 1/|/?| 2 <C Xq <C 1. In the first 
interval 

n xo , oo 



Pi(/3,x ) 



(3 J 77(1 

v 



X) 



exp(—(3r)z) 



\ sinh 2 



dz. 



(C.3) 
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Integrating over z by part and then over x we have 
„ .„ . 2x 8 °f dz , ^ , „. . „ , / coshz 1 \ 

ft(A *0 - ^ + W J ? (i - (/V* + i) «p(-/JV3S)) (-^ - ? J . 

(C.4) 

In the integral which contains the term f3z*Jx~Q the interval z C 1 contributes so 
that 

°° d 1 
j -^f3zy/x^exp(-f3zy/x^) y-j^J = (C5) 

o 

In the remaining integral we split the interval of integration into two: 1) < z < 
Zq-, 2)z < z < oo (1/(/3^/xq) -2 1), then 

dz ( cosh ,2 1 \ 

(1 - exp(-^ v /x )j 



z 1 \ sinh 3 z z 3 

o v 

zo 



f dz . . n , .. / z \ r dz I coshz 1 

J F (1 - exp(-^) (-- + / - ? 

20 x 

20 



(C.6) 



So we have 

Im P^o) - || - ^ (C.7) 

In the second interval over x (1 > x > xo) the interval zCl contributes as well 
as in the first term of P(/3) Eq.( |B.l| ) which we include here 

^0,^ = ^-1^ + 0(1) (C.8) 

Adding fltT7|) and (jU^) we get for > 1 

ImPW ^-W (a9) 
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Figure captions 



Fig.l The pair creation spectral probability p , x = — in gold in terms 

ax uj 

of the exact total Bethe-Heitler probability taken with the Coulomb cor- 
rections (see Eq. fl2.24j )). 

— Curve BH is the Bethe-Heitler spectral probability (see Eq.2.23); 

— curve Tl is the total contribution (the sum of the logarithmic ap- 
proximation dWp/de fl2.10|) and the first correction to the spectral 
probability dW^/de Q2.15| )) for the photon energy uj = 2.5 TeV; 

— curve 2 is the is the logarithmic approximation dW^/de ( [2.10| ), curve 
c2 is the first correction to the spectral probability dW^/de Q2.15Q ) 
and curve T2 is the sum of the previous contributions for the photon 
energy uj = 25 TeV; 

— curves 3, c3, T3 are the same for the photon energy uj = 250 TeV; 

— curves 4, c4, T4 are the same for the photon energy uj = 2500 TeV; 

Fig. 2 The relative energy losses of electron per unit time in terms of the 

Bethe-Heitler radiation length L® ad : -L° rad in gold vs the initial energy 

of electron (curve 1) and the total pair creation probability per unit time 
Wp (see Eq.( p.25| ))in terms of the Bethe-Heitler total probability of pair 
creation W^ H (see Eq.( [2.24| )) in gold vs the initial energy of photon (curve 
2). 

" o mi , • r ,. • dW dW UJ . 

b ig.3 the spectral intensity of radiation uj— — = x— — , x — — m gold m 

duj dx f 
terms of 3L ra d taken with the Coulomb corrections (see Eq.( |3.9|) ). 

— Curve BH is the Bethe-Heitler spectral intensity (see Eq.3.6); 

— curve 1 is the is the logarithmic approximation ujdW c /duj Eq.(2.28) of 
||, curve cl is the first correction to the spectral intensity udWi/duj 
Eq.(2.33) of [§] and curve Tl is the sum of the previous contributions 
for the electron energy e = 25 GeV; 

— curve 2 is the is the logarithmic approximation ujdW c /duj Eq.(2.28) of 
[§], curve c2 is the first correction to the spectral intensity udWi/duj 
Eq.(2.33) of [[J and curve T2 is the sum of the previous contributions 
for the electron energy e = 250 GeV; 

— curves 3, c3, T3 are the same for the electron energy e = 2.5 TeV; 

— curves 4, c4, T4 are the same for the electron energy e = 25 TeV; 
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